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Abstract: Collider experiments are one of the most promising ways to constrain Dark Mat¬ 
ter (DM) interactions. For several types of DM-Standard Model couplings, a meaningful 
interpretation of the results requires to go beyond effective field theory, considering simpli¬ 
fied models with light mediators. This is especially important in the case of loop-mediated 
interactions. In this paper we perform the first simplified model study of the magnetic dipole 
interacting DM, by including the one-loop momentum-dependent form factors that mediate 
the coupling - given by the Dark Penguin - in collider processes. We compute bounds from 
the monojet, monophoton, and diphoton searches at the 8 and 14 TeV LHC, and compare the 
results to those of direct and indirect detection experiments. Future searches at the 100 TeV 
hadron collider and at the ILC are also addressed. We find that the optimal search strategy 
requires loose cuts on the missing transverse energy, to capture the enhancement of the form 
factors near the threshold for on-shell production of the mediators. We consider both minimal 
models and models where an additional state beyond the DM is accessible. In the latter case, 
under the assumption of anarchic flavor structure in the dark sector, the LHC monophoton 
and diphoton searches will be able to set much stronger bounds than in the minimal scenario. 
A determination of the mass of the heavier dark fermion might be feasible using the Mt 2 
variable. In addition, if the Dark Penguin flavor structure is almost aligned with that of the 
DM mass, a displaced signal from the decay of the heavier dark fermion into the DM and 
photon can be observed. This allows us to set constraints on the mixings and couplings of 
the model from an existing search for non-pointing photons. 
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1 Introduction 

The existence of Dark Matter (DM) is firmly established by a large number of astrophysical 
and cosmological observations. Despite the fact that it contributes a large component of 
the energy density of the universe, however, its precise properties remain almost completely 
mysterious. The common belief is that most of the DM is in the form of a stable particle, which 
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is neutral or charged very weakly under the electric force, but interacts at least gravitationally 
with baryons. If such a particle carries non-gravitational interactions with the Standard Model 
(SM), the production of DM particles at high energy collider experiments offers one of the 
most promising opportunities to identify the nature of the DM interactions. Thanks to the 
capability to produce DM particles in a wide mass range and to the obvious independence 
from astrophysical uncertainties, collider searches provide complementary results to the direct 
and indirect detection experiments [ 1 - 6 ]. 

When setting collider constraints on DM interactions, the most straightforward way of 
parameterizing the DM-SM coupling is through an effective field theory (EFT) description. 
The non-observation of events with significant missing energy in excess of the SM background 
is then translated into upper bounds on the coefficients of the effective operators that couple 
the DM to the SM fields [7-21]. This simple method is independent of the details of the 
ultraviolet (UV) completion of the model. However, it is based on the assumption that the 
EFT gives a valid description of the collider process, i.e. that the mediators can be integrated 
out at the energy scale of the collision. Unfortunately, the sensitivity of the current LHC 
searches does not correspond to heavy enough mediators for many of the EFT couplings 
[22, 23], and unitarity usually sets stronger bounds than the collider search [24, 25]. Thus 
in many instances, to extract meaningful information from collider searches it is necessary to 
consider perturbative models with light mediators. This ‘simplihed model’ approach has been 
applied to several DM-SM couplings whose UV completions feature the tree-level exchange 
of mediators [26-39]. 

The inadequacy of the EFT approach is manifestly even more dramatic in the case of loop- 
mediated couplings. Among the loop-induced DM-SM interactions, an especially important 
status is held by the dipole and Rayleigh operators [40, 41], which play important roles in 
DM model building. Collider constraints on these couplings within the EFT approach have 
been explored in several previous studies [42-46], but the resulting bounds translate into very 
weak constraints on the mediator masses, thus calling for a simplihed model description. 

In this work we perform the hrst systematic collider study of a loop-induced DM-SM 
coupling, by considering a perturbative UV completion with light mediators and including 
the momentum-dependent form factors in the description of collider processes.^ In our study 
we focus on a generic UV completion of the magnetic and electric dipole operators, the Dark 
Penguin. We emphasize, however, that the method introduced in this work can be applied 
to any other loop-induced DM-SM coupling. 

If the DM is a singlet under the SM gauge symmetry, the magnetic and electric dipole 
operators fJ-M generated by the dark penguin diagrams, give the 

lowest order interactions of the DM with the SM gauge helds. As a consequence, these 
operators play important roles in the possible explanation of various 7 -ray excesses observed 
at indirect detection experiments [49-55]. Moreover, being chirality-hipping, these operators 

first brief discussion of the loop form factors for dipole-interacting DM production at colliders was 
presented in Ref. [47] . See also Ref. [48] , where the loop contribution of the dark sector to dilepton production 
was considered, although in a different model. 
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can possess a non-trivial flavor structure when more than one species of dark fermions Xi 
is present. In this case the photon dipole couplings in the mass basis can connect the light 
and excited DM state, and give interesting inelastic scattering signals at direct detection 
experiments [56-61]. The typical momentum exchange in the direct and indirect detection 
experiments is very small, therefore it is reasonable to use the EFT description even when 
the mediators are light with respect to collider energies. 

In order to allow a comparison to the results of direct and indirect detection experiments, 
the collider searches need to provide bounds on both the DM-mediator coupling and the 
mediator mass. The search for mediator decays in various UV completions of the dipole and 
Raleigh operators has been studied in [62]. However, since the decay of the mediators only 
depends on their branching ratios, the information about the DM-mediator coupling is lost 
in these searches. Thus the study of the dark penguin process would be crucial even if the 
mediators were directly discovered first. 

The dark penguin serves as a good example of a simplified model that can be constrained 
by different collider searches: depending on the ‘dark flavor’ structure - the flavor structure 
of dark fermions, - it can provide signals in the monojet, monophoton, diphoton, and even 
non-pointing photon searches. Besides allowing us to set meaningful constraints, the inclusion 
of light mediators helps us to identify the optimal cuts to be used in the collider searches, 
which are different from those commonly employed in the study of the EFT couplings. When 
the mediator mass is much smaller than the typical momentum exchange in the dipole, the 
missing transverse energy (MET) distribution of the DM signal is much softer than the 
one obtained assuming EFT couplings. Therefore the intuition of setting harder cuts to 
increase the signal excess no longer applies: on the contrary, the best strategy is to keep 
the cuts as low as possible, to include the enhancement of the form factors corresponding 
to the mediators being produced on-shell. It follows that in the search for the dark penguin 
with light mediators, lowering the background is more important than increasing the collider 
energy. Looking ahead towards future experiments, this implies that the International Linear 
Collider (ILC) would have the capability to set much stronger constraints than a very high 
energy (100 TeV) hadron collider, as will be shown in our analysis. 

The organization of the paper is as follows. In Sec. 2 we present the simplified model 
used in the analysis, and address the possibility of having displaced photon signals from an 
aligned flavor structure. We begin Sec. 3 by explaining our method for including the loop- 
mediated dipole couplings in the simulation of collider processes. Then we discuss in detail 
the missing energy searches used to set constraints on the simplified model, paying particular 
attention to the estimation of the systematic uncertainties in the projection to the 14 TeV 
LHC. In Sec. 4 we present the constraints from LEP, the 8 and 14 TeV LHC, and the future 
100 TeV collider and ILC on the dark penguin parameter space, assuming a single flavor 
of dark fermions. Furthermore, we point out that the sensitivity to light mediators can be 
improved by choosing MET cuts weaker than those used in the search for effective couplings. 
To give an example of the case with more than one dark flavors, we analyze the model with 
two dark fermions in Sec. 5, setting bounds from LEP and the monophoton and diphoton 
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searches at the LHC. For the diphoton channel, the possible application of the Mt 2 variable 
to determine the mass of the heavier dark fermion is discussed. In Sec. 6 we consider the 
dark penguin with a flavor structure almost aligned to the dark fermion masses, in which 
case displaced photon signals can be observed at the LHC. We use the 8 TeV ATLAS search 
for non-pointing photons to compute the bound on the dark mixing angle and coupling. In 
Sec. 7 the collider constraints on the dark penguin are compared to the current results from 
direct and indirect DM detection experiments, by showing the reach of the different searches 
on the magnetic dipole moment and the annihilation cross section of XX 77 ) respectively. 
We conclude by summarizing our result in Sec. 8. Finally, App. A contains general formulas 
for the dark penguin form factors, whereas App. B collects the basic statistics we used for 
setting limits. 

It is worth pointing out that Sec. 3 is somewhat technical, therefore the reader mainly 
interested in the results of our work might prefer, after having become familiar with the dark 
penguin in Sec. 2 , to move directly to Sec. 4. 

2 A simplified model of dark penguin 

Here we describe a simple UV completion of the magnetic and electric dipole operators, which 
will be employed throughout the paper. The model is similar to the one discussed in [49, 62] 
and contains dark Dirac fermions Xi with flavor index i = 1,..., Ny., a fermion mediator 'll;, and 
a scalar mediator (j). Both mediators carry some hypercharge Y, and the Xi s-re SM singlets. 
A specific assignment of dark charge which stabilizes the DM particle may affect the decay 
of mediators, but not the dark penguin process we are interested in. The Lagrangian in the 
mass basis is written as 

L D Xi{0 - mi)xi + ilP'ip -Mfill'll; + \(l>f + [Xf'ii; Hr Xi </> + Xf'ii’PL Xi4> + k-C-) , 

( 2 . 1 ) 

where = d^ — iYg'B^. Several simplifications will be made on the Lagrangian: we assume 
the fermion and scalar mediators to have the same mass, Mf = Mg = and also assume 
the Yukawa-type couplings Xi = X[ = Xf to be real. This in particular implies that no 
electric dipole moment operator is generated from Eq. (2.1). Depending on the details of 
a specific model, one of the mediators can decay into SM particles, while the mediator that 
carries the dark charge stabilizing the DM can decay into the DM and SM particles. A general 
analysis of the phenomenology of the mediators can be found in [62] . Their results show that 
the discovery prospects at the LHC depend strongly on the charge assignment. The most 
challenging scenario corresponds to SU ( 2 )£,-singlet mediators with Y = — 1 , and il; mixed with 
the right-handed tau lepton. Even at 14 TeV with 300 fb“^, the LHC will have no sensitivity 
to this model [62]. However, different hypercharge assignments can change dramatically the 
decays of the mediators, leading to significantly better prospects. In our analysis of the 
dark penguin we wish to be independent from the details of the model building, therefore we 

^Notice, however, that in App. A we present general results for the dark penguin with Mf ^ Ms. 
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Figure 1. Diagrams of the dark penguin. The corresponding amplitude is given in Eq. (2.2). 

treat V and M as free parameters. Precision electroweak measurements do not set relevant 
constraints on mediators carrying only hypercharge. Even if the mediators are doublets under 
SU{2)l, as considered for example in Ref. [49], no contribution to the S and T parameters 
arises at one loop. 

The model in Eq. (2.1) generates a magnetic dipole operator through the dark penguin 
process in Fig. 1. The amplitude for —)■ XiXj is written into a gauge invariant form 

^-^penguin = ^ QyF^] v{pj) , (2.2) 

In this paper, we use various collider searches to set an upper bound on^ 

XVyN. (2.3) 

where the factor AjAj encodes the dark flavor structure. The factor N is the multiplicity of 
the mediators {N > 1 can arise, for example, if the mediators transform non-trivially under 
SU (2) l). The form factors Fq^fj as functions of the masses (m^., m ^., Mf, Mg) and momentum 
(defined in Fig. 1) are given in App. A. It is important to note that the model discussed 
here can also generate other loops, such as those for the Rayleigh operator, . 

However, the contributions of this operator to the collider processes considered in this paper 
carry either extra gauge couplings or phase space suppressions, and thus only give sub-leading 
effects, except in a few cases on which we will comment in what follows. We leave the detailed 
study of the Rayleigh operators for future work. 

In presence of a non-trivial flavor structure of the dark fermions Xi^ important dif¬ 
ference between collider searches and (in) direct detection experiments is that the collider 
processes can generically involve more than one dark flavor. Indeed, as we will show, the 
current and upcoming missing energy searches at colliders are strikingly more sensitive to the 
case where more than one dark flavors are within kinematic reach. In this case, some assump¬ 
tions on the flavor structure are necessary in order to compare results between the different 

^Our analysis is insensitive to the sign of the hypercharge of the mediators. For simplicity, from now on we 
denote with Y the absolute value of the hypercharge. 
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experiments. In this paper, we assume the Aj couplings have a totally anarchic structure, 
with no unnatural hierarchies between different flavors. The assumption of anarchic struc¬ 
ture permits a direct comparison of the collider bonnds among themselves and to the direct 
and indirect detection experiments, which are typically sensitive only to the couplings of the 
lightest dark fermion (the DM particle). To be more precise, in the following we assume 
Aj = A (1 -|- dj) with |(5j| <C 1. For models with different flavor structures, one can rescale onr 
bounds from the different searches to set proper constraints. 

According to naive dimensional analysis (NDA), the perturbative bound is A\/iV < dvr. 
When showing our results, we allow X^/YN to be as large as Svr, which corresponds to the 
NDA perturbative bound for Y = 4 (corresponding to g'Y ~ 1). In a different UV completion 
of the dipole interaction, such as for example a model where the scalar mediator is replaced 
with a gauge boson, the quantitative result for the dark penguin would be modified. However, 
the fact that gauge invariance constrains the structure of the amplitude forbids a qualitative 
alteration of our analysis. 


2.1 Displaced signals and the aligned flavor structure 

Another interesting scenario to explore is the case where the flavor structures of the DM mass 
matrix and the dark penguin are nearly aligned. This means that the fermion mass matrix 
nfiijXiXj is almost proportional to the dark penguin, whose flavor structure is determined by 
Xi Xj 

mij oc XiXj{l + Cij) + rh6ij, \eij\ <C 1, (2.4) 

where the term proportional to the identity gives the lightest dark fermion a mass ~ m. When 
rotating the dark fermion fields into the mass basis, the dark penguin between different mass 
eigenstates carries an extra 0(e) snppression that can result into a displaced decay of a heavy 
dark fermion Xh to the light DM xi and a photon, with a width (assuming M 3> 3> rn^.^) 


F(xh xn) 


e^X'^Y‘^N‘^m\ „ 

X e^. 


87r(327r2)2M2 


(2.5) 


Contrarily to the case where the decay is displaced because of phase space suppression [63], 
here the photon is hard, therefore the process is accessible in LHC searches. 

The aligned flavor structure can be generated if the dark sector has a single flavor breaking 
spurion, and both the DM masses and dipole interactions are generated by the same loop-level 
mediation. Since the gange coupling is flavor blind, the flavor structures of the two operators 
are identical 


r\j 


A^ Vij eQ^ 
16 'k‘^M 






( 2 . 6 ) 


implying that no heavy-light dipole coupling is present in the mass basis for dark fermions. 
However, if there exists a small chiral symmetry breaking in the infrared, or an extra flavor 
symmetry breaking from an even higher order mediation is present, the misalignment between 
the newly generated mass and the dark penguin leads to a small heavy-light dipole coupling 
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and thus a suppressed decay rate as in Eq. (2.5). With the assumption of anarchic couplings 
Aj, the search of non-pointing photons-|-MET can set interesting bounds on the (A,e) plane. 
The result is discussed in Sec. 6. 

3 DM searches at colliders 

In this section we describe the collider analyses used in this paper to set bounds on the DM 
parameter space. We begin by outlining our procedure for taking into account the full loop 
dark penguin form factors in the collider simulations. Then we move on to describe in detail 
the relevant searches, including the monojet, monophoton and diphoton final states at the 8 
and 14 TeV LHC. Then we turn to the 8 TeV search for non-pointing photons. Einally, the 
projections to a 100 TeV collider and to the ILC are studied for some of the searches. 

In our analysis of 8 TeV searches, we reproduce the shape of each SM background using 
our MonteCarlo (MC) simulations, and compute the additional global rescaling factor needed 
to obtain exact agreement with the distributions reported in the experimental papers. We 
then apply these rescaling factors in the 14 TeV projections. In the 14 TeV projection of 
monojet and monophoton searches, we estimate the improvement of systematic uncertainties 
by separating them into two parts. For uncertainties that relate to the normalization of 
the SM background, we follow the data-driven analysis at 8 TeV by simulating the control 
region sample at 14 TeV. For other uncertainties, including those on the PDFs and the 
acceptance of the detector, we assume an improvement proportional to the square-root of the 
luminosity. We do not attempt to simulate the QCD background at 14 TeV in this work. In 
the 8 TeV search, the pure QCD processes contribute less than ~ 1% (10%) of the monojet 
(monophoton) background, and we make the reasonable assumption that the background can 
be kept subdominant at 14 TeV by setting a harder MET (photon pt) cut. On the other 
hand, the QCD background does play an important role in the 8 TeV diphoton-l-MET search. 
We include additional jet and lepton vetoes to suppress it at 14 TeV. 

We generated both the signals and backgrounds using MadGraph5 [64] and showered 
the parton level events using Pythia 6 [65]. We used PCS 4 for the detector simulation and 
cross-checked the resnlts using Delphes 3 [66]. To compute the QCD iL-factors for some SM 
backgrounds, we nsed MCFM [67] and VBFNLO [68]. 

3.1 Including loop form factors in a collider process 

To properly describe loop-mediated processes at colliders, we simulate the DM production 
using EFT operators and reweight the events by employing the expressions of the form factors 
given in App. A. As described in Eq. (2.2), the amplitude of the dark penguin contains three 
distinct Lorentz structures 

iFq{q'^) - m^,)q>^Bi,. (3.1) 

The coefficient of the last operator ~ vanishes in the Ny. = 1 case, since rriy.. = rn^y For 
Ny- > 1 the contribution of this operator to the amplitude for pair production of DM particles 
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// —>■ XiXj is proportional to Vf{pj, + pj)uf and is thus strongly suppressed by the light SM 
fermion masses, after using the equations of motion. Therefore we neglect the operator ~ 
altogether and only consider the and operators. 

For a given process, we first generate the events in MadGraphS using a linear combination 
of the two relevant effective operators, then we reweight the events using the ratio of the 
partonic matrix element squared computed retaining the form factors to the value obtained 
using the EFT. 

To give a concrete example, we consider the Ny. = 2 monophoton analysis. We first 
simulate the signal process in Fig. 13 using two operators implemented through FeynRules 
[69] 

Oa ~ ^Xh crf"’''xiB,,u + h.c., Ov ~ Xhl^XlB^, + h.c.. (3.2) 

After performing the detector simulation, we apply the cuts for the monophoton search and 
obtain a list of signal events. For each event, we obtain the value of q^ from the corresponding 
parton-level four momenta and use it to compute the form factors Fa,q{q‘^). We then reweight 
the event using the expression of the matrix element squared, schematically 

^<7dark penguin dfTEFT o\ 

"" \Mo.+Moyf ■ 

where ■M.FaqiQ'^) the amplitudes corresponding to the first and second term in Eq. (3.1), 
respectively, whereas Aio^v the amplitudes corresponding to the effective operators. 

The procedure described above fully accounts for the interference between the two rel¬ 
evant Lorentz structures, and was applied in the study of the iV^ = 2 scenario, where both 
monophoton and diphoton signals are given by 2 —2 scatterings followed by the decay of Xh, 
as shown in Eigs. 13 and 16. On the other hand, for = 1 the main constraint comes from 
the monojet process in Eig. 10, which is genuinely a 2 —)• 3 scattering. Eor the sake of sim¬ 
plicity, in this case we neglect the interference term in Eq. (3.3), which carries an additional 
~ niy./^/q^ suppression due to the different chirality structure between the two operators. In 
the Ny. = 2 case, the inclusion of the interference term gives at 14 TeV an increase of the cross 
section of less than 40%, which translates into a correction of less than 10% to the constraint 
on AVT N. Such a deviation is acceptable when compared to the possible uncertainties in 
our projection. Thus we simulate the monojet events using the incoherent sum of the two 
operators in Eq. (3.2) to obtain the dark penguin result. 

To demonstrate the effect of the form factors, we plot in Fig. 2 the distribution 
of the monojet events in both the dark penguin and EFT cases. Notice that due to the 
kinematics of the process, the MET is given by \/(p multiplied by some angular factors. The 
distribution of the dark penguin is remarkably different from the one of the EFT: the dark 
penguin exhibits an enhancement around ^/q^ ~ 2M, corresponding to the threshold for the 

"^The sum of the second and third term in Eq. (3.1) actually corresponds to the gauge invariant dimension-6 
operator /h? -|-h.c. (see also App. A). To simplify the event generation we use Ov instead, and 

include the factor from the derivatives in the reweighting of the events. 








Figure 2. The normalized distribution of the DM signal in the monojet search. The blue his¬ 
tograms correspond to the simulation of the effective couplings xB^.v/^ (left) and Xl^BjIs? 

(right), while the pink (orange) histograms correspond to the full dark penguin form factors with me¬ 
diator mass M = 300 (500) GeV. The form factors peak at \/^ ~ 2M. The cuts applied are described 
in Sec. 3.2 (the cut on the MET is set to 550 GeV). 

production of on-shell mediators [47], and it is softer than the EFT at high energy. As a 
consequence, contrarily to the strategy used for effective DM couplings, where a strong cut 
on the MET is generically preferred, in the search for the dark penguin with light mediators 
a softer cut is favored. We will return to this point when discussing the monojet result, see 
Fig. 12. The preference for soft MET cuts applies also to the monophoton and diphoton 
searches. 

3.2 LHC monojet 

For the monojet channel, we follow and extend the CMS analysis in [70]. The event selection 
requires one jet with prUi) > HO GeV and \rj{ji)\ < 2.4. A second jet with PTij 2 ) > 30 
GeV and \r]{j 2 )\ < 4.5 is allowed, as long as the A(j){ji,j 2 ) < 2.5. Events containing a third 
jet satisfying priJs) > 30 GeV and \r]{j 3 )\ < 4.5 are vetoed. To reduce the backgrounds from 
Z and W production, events containing electrons with pT > 10 GeV and jr/j < 2.5, muons 
with pt > 10 GeV and \r]\ < 2.1, or hadronic taus with pT > 20 GeV and \t]\ < 2.3 are 
also vetoed. The counting experiments in [70] are performed in 7 signal regions, with MET 
cut from 250 to 550 GeV with a step of 50 GeV. In the left panel of Fig. 3 we show the 
comparison of our simulated backgrounds with those reported by CMS. For each of the three 
main backgrounds Z-|-jets, IE-|-jets and tt, we use MCFM to obtain the QCD A'-factor and 
apply an additional rescaling factor to match the CMS result. The rescaling factors are 0.84, 
0.95 and 0.60, respectively. The rescaling factor for Z-|-jets is also applied to the DM signal 
events. 

The CMS search uses a data driven analysis based on a /i-|-jets control sample to determine 
the normalization of the jZ and jW backgrounds. The corresponding events in the control 
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Figure 3. Left panel: the dominant SM backgrounds at 8 TeV from the CMS study in [70] (shaded 
regions) and from our MC simulation (dashed lines). Right panel: distribution of the dominant 
backgrounds at 14 TeV. Here we take a large DM coupling, \\/YN = 30, for the visualization of the 
signal distribution. 




Figure 4. Event distributions of Z’(/i+/r )+jets and lF(^i^)+jets in the control region. 

sample are ^~)+]eis and VF(/«^)+jets, respectively. The cuts for the control sample 
are the same as for the monojet search, except the lepton vetoes are not applied. The 
Z{^'^ ^~)+ieis sample requires two muons with pT > 20 GeV and \t]\ < 2 . 1 , with at least 
one of the muons passing isolation requirements, and the invariant mass of the muon pair 
between 60 and 120 GeV. Similarly, the lT(/izz)+jets sample requires one isolated muon with 
Pt > 20 GeV and |r/| < 2.1, and the transverse mass of the muon plus neutrino system in the 
range 50 GeV < Mt < 100 GeV. Fig. 4 shows the comparison between our MC simulation of 
the control regions and the CMS results. For Z{p'^p,~)+iets, the MET is defined as the sum 
of the muon transverse momenta, while for lT(/izz)+jets the MET is given by the neutrino 
transverse momentum. The good agreement with the CMS results allows us to simulate the 
data driven analysis in the 14 TeV study. 

Eor the 14 TeV projection we follow the same cuts in the CMS 8 TeV analysis, apart from 
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Figure 5. Monojet and monophoton systematic uncertainties as functions of luminosity and MET 
cut, following the assumptions in Tables 1 and 2. The number shown on each contour is the total 
systematic uncertainty divided by the number of background events. 



CR statistics 

Other systematics 

Z+jets 

^^CRzi 

3.9 £-1/2 

IT+jets 

^^CRwi 

4.6 £-1/2 


Table 1. Summary of the contributions (in %) to the uncertainty used for the 14 TeV monojet study, 
following the analysis in [70]. The dominant uncertainty comes from the limited number of control 
sample events. £ = £/(20 fb“^), and we assume the other sources of uncertainty will be improved 
with the increase of data. 

varying the MET cut from 550 to 2250 GeV with a step of 100 GeV. The dominant systematic 
uncertainty for each MET cut choice is obtained from the control sample simulation. For the 
remaining uncertainties, we take the values quoted in the GMS analysis for a MET cut of 
550 GeV, and assume they will decrease with the square root of luminosity, see Table 1. 
Since at 8 TeV these additional uncertainties do not depend strongly on the MET cut, it 
is a reasonable guess to use their value also at 14 TeV. The right panel of Fig. 3 shows the 
projected monojet backgrounds studied in this work, together with the dark penguin signal. 
The total systematic uncertainty is shown in Fig. 5 as a function of luminosity and MET 
cut. To provide some figure of merit, with 3 ab“^ of data and MET cut of 550 GeV the 
background is ~ 10® events, which requires ~ 10^ DM signal events for a few-cr excess, given 
an uncertainty of order 1%. The corresponding \^/YN is ~ 20 when the mediator mass is 
M ~ 500 GeV. 
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Figure 6 . Left panel: comparison of the 8 TeV monophoton backgrounds from our MC simulations 
and from the CMS analysis. Our results are shown as dashed histograms, stacked on top of the sum of 
the subleading backgrounds Z{ll)j, 7 j, 77 , QCD and beam halo) as given by CMS [71], in 

order to facilitate the comparison. Right panel: the projected 14 TeV background used in the analysis. 

3.3 LHC monophoton 

For the monophoton channel, we follow the 8 TeV CMS search in [71]. The event selection 
requires one photon with > 145 GeV and < 1.4442, and in addition > 140 GeV 
with > 2. Events containing electrons or muons with px > 10 GeV, jry] < 2.5, and 

> 0.5 are vetoed, as well as events containing more than one reconstructed jet with 
px > 30 GeV, \r]\ < 2.5 and ARj^ > 0.5. The comparison of our simulation of the photon px 
distribution to the one obtained in the CMS analysis is given in Fig. 6 . For each background 
we apply a p^-dependent iC-factor obtained from the NLO MCFM calculation [67].® In addi¬ 
tion, we apply an overall rescaling in order to match the normalization provided by CMS: the 
rescaling factor is equal to 1.1 for Z{i?v)'j, 1.2 for W{ev) and 1.8 for W{lv)'y.^ We do not at¬ 
tempt to simulate the backgrounds given by the misidentification of leptons {W{pLiy), Z{ll)'y), 
jets ( 77 ), 77 , beam halo, and QCD, but these backgrounds are subdominant as is shown in 
the left plot of Fig. 6 . Only the three dominant backgrounds Z{vu)^, IF(ez^), are 

included in the 14 TeV analysis. 

In the CMS analysis [71] , the main source of systematic uncertainty is given by the higher- 
order QCD corrections to the Z{vv)^ and W{^v)^ backgrounds. Although the collaboration 
performs a data-driven analysis for each of the two backgrounds and compares the results to 
those obtained from MCFM, the limited number of control sample events at 8 TeV does not 

®The K-iactor is very large for W(lv)Xi K > 4.5 in the region considered > 145 GeV and growing with 
increasing p^. This is a consequence of the presence of a ‘radiation zero’ [72, 73] that strongly suppresses the 
LO amplitude. The NLO QCD corrections (in particular those associated with the qg partonic channel, which 
is not suppressed by the radiation zero) therefore contribute a large fraction of the total cross section, leading 
to the big Jf-factor [74]. 

®We do not apply to the DM signal the rescaling factor 1.1 obtained from the Z{vv)'y background, since it 
amounts to a negligible correction. 
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Figure 7. Control regions for the Z 7 (left) and (right) backgrounds, with the inverted angular 
cut < 2.9 for the Z(yv)^ events and the inverted lepton veto for the events. 

help to reduce the uncertainty. Therefore, in our 8 TeV analysis we made directly use of the 
uncertainties quoted by CMS. However, following the same strategy of the monojet case, we 
are going to assume that the data-driven analysis will play an important role in the future. 
With the increase of luminosity, the systematic uncertainty will be based on the statistics of 
the control sample. In this case, using too stringent kinematic cuts in the search can increase 
the size of the systematic uncertainty, and a proper choice of the cut is required to optimize 
the signal significance. 

To perform a data-driven analysis at 14 TeV, we compare our 8 TeV simulation of the two 
control region samples to the CMS results, and then apply the same analysis for the 14 TeV 
case. The Z{vv)^ control region is defined by the inverted angular cut A(j){^rp, 7 ) < 2.9, while 
the W{iiy)'y control region is dehned by the inverted lepton veto. In both cases we obtain 
a result in good agreement with CMS, see Fig. 77 This allows us to simulate the control 
region sample at 14 TeV and derive the relative systematic uncertainty shown in Fig. 5. As 
given in Table 2, the total uncertainty also contains various systematic uncertainties in the 
control region analysis for Z 7 and IT 7 , as well as the uncertainty on the probability for an 
electron to be misidentihed as photon for the IT(ez^) background. We take the values of 
these uncertainties from the 8 TeV analysis [71], and assume they will be improved with the 
increased luminosity as This is based on the fact that various data-driven analyses 

have been used to determine these uncertainties at 8 TeV. 

For the 14 TeV monophoton projection, we follow the same event selection used in the 
8 TeV CMS search, except for tighter cuts on the photon px and MET: both are required 
to be larger than 300 GeV, to suppress the QCD and beam halo backgrounds, which are not 
included in the 14 TeV estimate. As discussed in Sec. 3.1, choosing a tighter cut does not 
improve the sensitivity to the dark penguin. 

^Notice that we need to apply large rescaling factors for the W(lv)^ and W{ev) backgronnds in the Z 7 
control region, 4.2 and 3.0, respectively. However, these backgrounds are subdominant to Z{vi')y. 


- 13 - 






























Background 

Relative systematic uncertainty 

Z{vv)x 

\INcr,., + ( 0 . 10 £- V 2 iV ^^^^)2 + 

W{lv)x 

\jNcRtot + 

W{eiy) 

0 . 10 £- i /2 


Table 2. Summary of the contributions to the uncertainty used for the 14 TeV monophoton search. 
£ = £/(20 fb“^), and the number of control sample events depends on the photon pt cut. The relative 
uncertainties are taken from the 8 TeV CMS analysis [71]. 
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'Pjet 
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Ht > (GeV) 

400 

600 
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0 

^miss y (QgY) 

200 
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jet veto 
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no veto 

no veto 

PT > 50 GeV, \r]\ < 2.8 

lepton veto 

no veto 

no veto 

no veto 

PT > 25 GeV, \r]\ < 2.5 

Predicted 

1.01 + 0.36 

2.38 + 0.69 

1.59 + 0.58 

2.90 

Observed 

1 

5 

2 

- 


Table 3. Cuts, number of predicted events and number of observed events for the signal regions 
used in the 8 TeV ATLAS diphoton+MET analysis and in our 14 TeV projection. For the latter a 
luminosity of 300 fb“^ was assumed. 

3.4 LHC diphoton+MET: prompt 

If there is more than one flavor of Xi > I); ^ heavier Xh can decay to the lightest xi 
a photon: Xh —t xil- If Xh is pair produced at the collider, diphoton+MET can be important 
in probing the multi-flavor scenario. For this channel, we follow the diphoton search at 8 TeV 
of energy and 20.3 fb“^ of data performed by the ATLAS collaboration [75]. The search is 
aimed to constrain the parameter space of gauge mediated supersymmetry breaking models. 
The event selection applied in this analysis is as follows. At least two photons with > 75 
GeV and \r]^\ < 2.37 are required. Jets are reconstructed using the anti-fet algorithm with 
a radius parameter of 0.4. The jets are required to have pj, > 30 GeV and \r]j\ < 2.8. No 
vetoes are applied on the number of leptons and jets. Several new variables are introduced 
as follows. An angular separation variable, is defined as the minimum azimuthal angle 

between and the two selected photons. In presence of jets, a variable is introduced 

and defined as the minimum azimuthal angle between E™'®® and the two highest reconstructed 
jets, where the jets are required to have pj, > 75 GeV. The total visible energy, Hp, is defined 
as the scalar sum of the transverse momenta of photons, jets and leptons. 

There are three relevant signal regions defined in the first three columns of Table 3. 
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The region WPl is preferred for a larger mass splitting between and rn^^, because the 
cut Ht > 400 GeV can only be satisfied with high pT photons. In this signal region the 
SM background is smaller than in the other two regions, therefore a stricter bound can be 
achieved. For a lower mass splitting, the region MIS is preferred because it does not have any 
requirement on the Ht value. At 8 TeV, the region WP2 provides a weaker bound compared 
with the other regions, because of an upward fluctuation in the number of observed events. 

Note that none of the signal regions veto on the leptons and the jets. In the benchmark 
model used in the ATLAS analysis, the NLSP are produced in the decay chain of either gluino 
or chargino, and either jets or leptons are always present in the final states. At the partonic 
level, the final states of our model contain only a pair of photons and missing energy from 
the DM particles xi- Hence vetoing jets and leptons increases the sensitivity of the search to 
our model. 

The main backgrounds for these regions are IT 77 , Zjj and “QCD”. In this case, “QCD” 
is defined as the sum of multi-jet, 7 + jets and 77 + jets processes. When estimating the 
background, we multiplied the LO cross section of the IT 77 (^ 77 ) background obtained from 
MadGraphS with a iP-factor of 8.1 (1.8) obtained from VBFNLO [ 68 ]. This result still needs 
to be multiplied by a factor of 3.7 (1.1) to match the distribution reported by ATLAS.® The 
sum of the IT 77 and Zjj missing energy distributions in the WP2 signal region is shown in 
Fig. 8 . In this region the IT 77 background dominates Z^x by an order of magnitude. 

For the 14 TeV simulation, we use cuts inspired by the MIS region of the 8 TeV ATLAS 
analysis, see Table 3. The signal region is defined to have a slightly tighter missing tranverse 
energy than the 8 TeV analysis. Additionally, in order to increase the sensitivity of the 
search, jet and lepton vetoes are applied. In order to take into account the low pT jet veto 
correctly, we generated jet matched samples for both the signal and backgrounds. Gontrarily 
to the 8 TeV case, at 14 TeV Z 77 is the leading background, as a consequence of the lepton 
veto: the expected background for the 14 TeV signal region is of 2.64 events from Zjj and 
0.26 events from IF 77 , assuming 300 fb“^ of luminosity. In the 8 TeV ATLAS analysis, 
the Zxx background was estimated using MG, with an associated systematic uncertainty of 
50%. Similarly to the previous subsections, we assume that the systematics will improve 
with luminosity as The systematic error is thus estimated to be of 0.69 events at 

300 fb“^ and 2.2 events at 3 ab“^. We neglect the systematics for the subdominant IT 77 
background, since the number of events expected from IT 77 is smaller than the uncertainty 
on Z 77 . Finally we did not simulate the “QCD” backgrounds, but they are expected to be 
small due to the jet and lepton vetoes together with a harder MET cut. 

®In analogy with the IT 7 process, the very large N-factor for VL 77 can be understood as due to the presence 
of a radiation zero in the LO amplitude when the two photons are collinear [76]. We cannot explain the large 
rescaling factor of 3.7 needed to match the ATLAS IF 77 result, which was obtained by means of a data driven 
analysis. 
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Figure 8. Comparison between the ATLAS result and our MC simulation for the sum of VF 77 and 
Z 77 in the signal region WP2 (described in Table 3). 


3.5 LHC diphoton+MET: displaced 

For the search of the displaced photon signal plus missing energy, we follow the nonpoint¬ 
ing photon analysis in [77], performed by the ATLAS collaboration on about 20 fb“^ of 8 
TeV data. The full search also uses the delayed photon measurement, however, due to the 
complication of modeling the time of flight of the photon from the displaced vertex to the 
electromagnetic calorimeter (ECAL), we only focus on the measurement of Az-,, of nonpoint¬ 
ing photons (see Fig. 9). For DM signals given by the long-lived Xh —^ A/7 decay, /S.z^ can be 
related to the Xh decay length id in the lab frame: 


Az.y = 4 




rxh,T-r^,T^ 

,2 'T ,2 


1 - (^7,2)^ 




cos — cos((j(>^^ — (pj)cot sin 6 


Xh 


(3.4) 


where tt^z represent the transverse and longitudinal components of the unit vector r, re¬ 
spectively, as shown in Fig. 9. To obtain the distribution of the DM decay, we first 
simulate the prompt process, pp —>■ XhXhiXh XiliXh A/7 in MadGraphS, apply the 
cuts performed in the ATLAS analysis, and reweight the events using the dark penguin form 
factors. Then we calculate the proper lifetime of Xh and boost it to the lab frame using the 
momenta of each parton-level event. The angular information of the photon and Xh allow us 
to calculate Az^ in Eq. (3.4) as a function of the decay length. Using this, each simulated 
MC event contributes to the differential cross section in Az^ as 


<^<Tciispiaced _ _ l/^l —iiAz-i 

dAz^ ~ “ '^P^^P* 2 ® 


(3.5) 
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where the /j. characterizing the probability distribution dP/dAz^ of the decay is defined as 

rxh,T • 


— ^Xh^'^Xh 

h — 

PXh 


rxh,z - 


1 (^7,^) 
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'7 

-1 


(3.6) 


Summing the distributions derived from all the simulated events we obtain the differential 
cross section in Az-,,, see Fig. 9. 

The ATLAS search requires at least two loose photons with |r/| < 2.37 and Et > 50 GeV. 
At least one photon is required to be in the barrel region \r]\ < 1.37. To avoid collisions due to 
satellite bunches, both photons are required to have an arrival time at the ECAL t-f smaller 
than 4 ns, with zero defined as the expected time of arrival for a prompt photon from the 
primary vertex. We do not attempt to fully simulate which would require a more complex 
detector description, but rather we approximate with the time of flight of the Xh, requiring 
it to be smaller than 4 ns. In our estimation we do not include the detailed isolation cuts on 
the photon. We also neglect the effect of the displaced decay on the angular acceptance of 
the photons, simply imposing the requirements on \7]\ at the level of the prompt event. The 
signal region also requires > 75 GeV. Finally, to simplify the discussion we assume that 
every event has a reconstructed primary vertex in the geometrical center of the detector. 

For events where only one photon satisfies \r]\ < 1.37 (i.e. it is in the barrel calorimeter), 
this photon is used for the measurement of Az-f. For events where both photons are in the 
barrel, the photon with larger t-y is used. We approximate this timing condition by taking 
the photon emitted by the more boosted Xh, in which case the average decay is more delayed. 
In Fig. 9 the generated Az^ signal distribution is shown, on top of the expected background. 
The latter is taken from Fig. 4 of the ATLAS paper [77]. Because we are focusing on the 
non-pointing photon signals, to set constraints on the DM couplings we remove events with 
|A 2 :-),| < 30 mm. In our exploratory analysis we only consider the statistical uncertainty on 
the background, neglecting the effect of systematics. 


3.6 100 TeV collider and ILC 

To give an idea of how much higher energy colliders can improve the sensitivity to the dark 
penguin, we estimate the bound from the mono-jet search at a 100 TeV collider. Following 
the discussion in [78], we simulate the DM signal and backgrounds with ^/s = 100 TeV and 
impose the cuts pT > 2.5 TeV, \t]\ < 2.2 on the leading jet, together with $rp > 3 TeV. 
A second jet with pT > 100 GeV is allowed as long as it has \r]\ < 4.5 and the azimuthal 
separation from the leading jet is A4> < 2.2. Events with leptons (taus) are vetoed if the 
lepton satisfies \r]\ < 2.5 and px > 20 (40) GeV. We assume a 2% systematic uncertainty 
when calculating bounds with 3 ab“^ of data. 

To show the importance of lowering the SM background for the dark penguin search, 
we also estimate the monophoton bound from the ILG-500P and ILG-IOOOP scenarios in 
Sec. 4. Our analysis follows the one in [79] by assuming a 500 GeV (1 TeV) ILG with 250 
(500) fb“^ of data, and polarizations equal to P_ = -1-0.8 and = 0.5 for the electron 
and positron beams, respectively. The Z-related SM background e+e“ —)• Z{vi>)x can be 
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Figure 9. Left panel: the distribution of the non-pointing photon signals measured by ATLAS. 
The background reported by ATLAS (blue histogram) was obtained from a data driven analysis, using 
diphoton events with < 20 GeV. Also shown, stacked on top of the background (red histogram), 
is the signal distribution for {m^^, , M) = (300, 10, 300) GeV, XVNY = 6, and e = 10“^. Right 

panel: the geometry of the displaced signals. 

eliminated by cutting away the Z pole, but the background process involving a t-channel W 
cannot be reduced by simple kinematic cuts. The polarization of the ILC beams plays an 
important role in reducing the latter background, since the W only couples to the left-handed 
electron. The proposed search requires > 8 GeV, |cos0.y| < 0.995, and imposes a veto on 
events with photon energy 238 < < 245 (490 < E^ < 495) GeV at ILC-500P (-1000P) in 

order to suppress the Z[vv)^ background. Notice that in our simulation of the dark penguin 
monophoton signals at the ILC, we only considered the case where the photon is emitted by 
ISR. However, in principle the contribution from the Rayleigh operator is of the 

same order and should be consistently included. This could lead to quantitative changes in 
our results, although the qualitative features would remain the same. We leave the inclusion 
of the Rayleigh operators for future work. 

4 Collider phenomenology of the = 1 case 

In this section we present the collider reach on the scenario where only a single species of 
dark fermion x is accessible at colliders. In this case, signals in the monojet and monophoton 
plus missing energy channels can be produced through initial state radiation. In addition, 
when the DM mass is less than half of the Z mass, the measurement of the invisible Z decay 
at the Z-pole also sets a strong constraint. 

4.1 Bound from the Z invisible width 

If the DM is lighter than vnzl‘i‘^ the decay Z —)> XX is kinematically allowed and constrained 
by the measurement at LEP/SLD of the invisible decay width of the Z. For mz 'C M, the 
decay Z —)> xX is well described by an effective dipole interaction xcr^'^xB^iu (see Eq. (A.12)), 
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leading to the following approximate expression for the width 


r(z xx) ^ 



24576 ttS M2 


^1 + 8 



4m2 


(4.1) 


The uncertainty on r(Z —>■ invisible) is of 1.5 MeV [80], therefore the 95% bound reads 


BR(Z ^ yy) < 1 X 10-3 . 


(4.2) 


This bound is shown as an orange curve in the right panel of Fig. 11. 

4.2 LHC and future colliders 

At the LHC, the most promising search for the dark penguin with = 1 is monojet, in 
which a high pT jet is produced through the QCD ISR process, see Fig. 10. The details of 
the analysis were given in Sec. 3.2. The signal and background yields and the systematic 
uncertainty on the background all depend on the MET cut and luminosity. The projected 14 
TeV constraint is computed varying the MET cut from 550 to 2250 GeV with a step of 100 
GeV, and taking the strongest bound. The results are shown as solid blue curves in Eig. 11 
for = 50 GeV (left panel) and 10 GeV (right panel). The LHC reach is limited; the 14 
TeV run can only cover part of the perturbative parameter space for relatively large V ~ 4, 
in which case the perturbative bound is A\/VN < Svr. 

To compare the dark penguin to the EET, we also compute the bounds by parameterizing 
the DM-SM coupling with the effective dipole interaction. The resulting constraint, shown 
by the dashed blue lines in Eig. 11, is weaker than the dark penguin one, since the dipole 
EET neglects the enhancement of the form factor for -\/^ — 2M, as well as the sizable 
contribution from the term in Eq. (2.2) when is large (in the EET the q'^j^ piece 
would correspond to a dimension-6 operator, see the first term in Eq. (A.12)). Eurthermore, 
the optimal MET cut for the dark penguin signal is typically softer than for the EET coupling. 
In Eig. 12 we show a comparison of the best MET cut for the dark penguin and the effective 
dipole coupling, as a function of the mediator mass. When the mediators are light, the best 
MET cut for the dark penguin is lower than for the effective coupling. Notice that in both 
cases the optimal MET cut increases at larger luminosity, because in our data-driven analysis 
the systematic uncertainties are mainly determined by the statistics of the control samples. 

The LHG monophoton search gives a much weaker bound compared to monojet, and will 
not be discussed in detail. 

Given the limited sensitivity of the LHG to the = 1 scenario, it is important to 
estimate how the bounds could be improved at future colliders. The study of DM production 
at very high energy machines, such as a 100 TeV pp collider, mandates the full inclusion of 
the loop form factors, since the typical partonic energy is much larger than the masses of the 
mediators. Differently than in the search for heavy particles, however, the significance of the 
dark penguin signal will not be appreciably improved at a higher energy machine, unless the 
systematic uncertainty can be greatly reduced. In fact, when the mediators have a sub-TeV 
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Figure 10. Diagrams for the monojet process. The grey circle indicates the dark penguin. 




M (GeV) M (GeV) 


Figure 11. Upper bounds on the DM-mediators coupling for = 1, as a function of the mediator 
mass M. The dark matter mass is assumed to be = 50 GeV (left) and 10 GeV (right). Here we 
optimize the LHC bound by choosing the best MET cut for each mediator mass. The blue dashed 
curve shows the bound from the effective dipole interaction, in Eq. (A. 12). The monophoton 

constraint from the ILC500-P (ILCIOOO-P) assumes 250 (500) fb“^ of data with 500 GeV (1 TeV) center 
of mass energy and a polarization P_ = +0.8 and P+ = 0.5. The orange curve shows the constraint 
from the current invisible Z decay measurement. The perturbative bound on Xy/YN depends on the 
hypercharge of the mediator, X\/YN < Itta/V. 


mass, both the SM background and DM production at 100 TeV are mediated by effectively 
nearly-massless particles. This results in similar shapes of the MET distribution for the signal 
and background, and as long as the background is systematics-dominated, the ratio between 
the signal and background does not vary significantly when increasing the collider energy. 
The pink curves in Fig. 11 give an estimate of the 100 TeV reach (the search is described in 
detail in Sec. 3.6). The result is an improvement of the bound on A\/V N by only a factor ~ 2 
with respect to the projected 14 TeV LHC constraint. The sensitivity of the 100 TeV collider 
ameliorates if the MET cut can be kept very weak, below the TeV scale, as in this case the 
signal region includes the enhancement of the form factors for ~ 2M, thus increasing 
the signal yield. However, it is not clear whether such a low MET cut would be achievable, 
and if so, whether the systematic uncertainties would be increased beyond the 2% that we are 
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Figure 12. The MET cut of the projected 14 TeV search that optimizes the DM signal significance 
as a function of luminosity. As can be seen in the plots, dark penguins with light mediators prefer 
lower MET cuts. When the mediators are heavier the MET distribution becomes harder, similar to 
the EFT case, and a tighter cut is preferred. The optimal MET cut increases with luminosity, because 
in our data-driven analysis the systematic uncertainties are mainly controlled by the statistics of the 
control samples. 

using as benchmark in our estimate. More detailed studies are necessary to precisely assess 
the reach of a 100 TeV machine. 

In comparison to hadron colliders, the monophoton search at the ILC appears to be 
more promising for the dark penguin signals. As discussed in Sec. 3.6, thanks to the beam 
polarization and full reconstruction of the missing invariant mass, the SM background can 
be efficiently reduced. As is shown in Fig. 11, the ILC500-P with 500 GeV center of mass 
energy, 250 fb“^ of data and polarized beams can improve greatly the sensitivity to the dark 
penguin. 

5 Collider phenomenology of the = 2 case with anarchic DM conplings 

Differently from the direct and indirect detection experiments, colliders are capable of probing 
extended dark sectors, through the production of the additional states that accompany the 
DM. In this section we study the monophoton, diphoton, and rare Z decay constraints on 
the = 2 dark penguin model, where one heavy dark fermion Xh is present in addition to 
the light DM particle xi- this scenario, large photon signals can be generated through the 
pair production of dark fermions, followed by the decay Xh —^ A/7- Under the assumption of 
anarchic couplings between the dark fermions and the mediators, we can compare the bounds 
from different searches and provide complementary information to the direct and indirect 
detection experiments. In the diphoton search, the kinematic variable Mt 2 can be used to 
determine the mass of Xh- 

The constraint from the invisible Z decay width applies also in the N-^ = 2 case, as long 
as m^^ < mzl‘2, as we assume. In addition, when also < mzl‘i^ the search at LEPI 
for events containing two photons plus missing energy sets a much stronger constraint on the 
DM couplings. Therefore we focus first on the case > 771^/2, and discuss separately the 
more constrained possibility of a lighter Xh, in Sec. 5.3. 
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Figure 13. Dominant process for the monophoton signal. The grey circle indicates the dark penguin. 

5.1 LHC monophoton 

The monophoton signal arises through the production of XhXU followed by the decay Xh 
Xn, see Fig. 13. The branching ratio for the decay is either ~ ~ 0.8 or unity, depending on 

whether the xi^ channel is kinematically open or closed. Therefore the cross section is much 
larger than for the ISR process discussed in the case Ny. = 1, whose amplitude is comparatively 
suppressed by the three-body phase space and by the additional electromagnetic coupling. In 
this section we calculate the bound on \y/YN from the 8 and 14 TeV monophoton searches, 
and study the dependence of the bound on the mass splitting between Xh and xi ■ The result 
is shown in the left panel of Fig. 14. For light mediators, the slope of the bound on X\/YN 
decreases. This can be understood by noticing that the requirement of a minimum photon 
transverse momentum p7p imposes while the dark penguin peaks at ~ 2M. 

Therefore, for M < p)p a, significant fraction of the signal does not pass the event selection. 
The effect is more evident in the 14 TeV case, where the cut on the photon pT is stronger 
(300 GeV). 

To gain some insight about how the mass difference between Xh and xi affects the 
monophoton signal, we fix the DM mass to = 10 GeV and show the photon pT dis¬ 
tribution for various choices of rriy-^ in the right plot of Fig. 14. The coupling A\/V-V is fixed 
to the value corresponding approximately to the reach of the 14 TeV LHG. We see that the 
signal yield decreases significantly only for < 20 GeV. Thus the monophoton search 
has sensitivity even for fairly degenerate dark fermions, with mass splitting — my.j much 
smaller than the required photon px, thanks to the boost of Xh from the production process. 

In Fig. 15 we show a comparison of the bound obtained using the dark penguin and the 
one computed using the effective dipole coupling Xh^^^^XiB^v + b.c. (see Eq. (A.12)). For 
each mediator mass we vary the photon px cut between 300 and 1000 GeV in 100 GeV steps, 
and choose the value that gives the strongest constraint (therefore the dark penguin bound 
shown in Fig. 15 is slightly better than the one in Fig. 14, which was obtained with a fixed cut 
> 300 GeV for all mediator masses). Similarly to the monojet case, the bounds computed 
by taking into account the full dark penguin are stronger than the EFT results, due to the 
form factor enhancement at \/^ ~ 2M and the contribution of the term in Eq. (2.2) 
for large (f'. 
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Figure 14. Left panel: LHC monophoton bounds for = 2. Right panel: dependence of the photon 
Pt distribution at 14 TeV on the mass of Xh- To generate the distribution, the cuts on the photon 
transverse momentum and MET were relaxed to > 100 GeV. The luminosity is fixed to 300 

fb“^. The vertical dashed line indicates the actual cut applied in the analysis, p^ > 300 GeV. The 
perturbative bound on Xx'YN depends on the hypercharge of the mediator, Xx'YN < 47r\/Y. 



Figure 15. A comparison between the monophoton bounds computed using the dark penguin ampli¬ 
tude (blue) and the effective dipole interaction Xh<y^'^XiB + h.c. (orange). For each mediator mass, 
the bounds are optimized by choosing the best photon pf cut between 300 and 1000 GeV. 

5.2 LHC diphoton-l-MET 

If the Xh is light enough to be pair produced with sizable rate, the diphoton-l-MET signal can 
be generated through the process shown in Fig. 16. Differently from the monophoton case, 
there can be a partial cancellation of the transverse momenta of the two xi particles, leading 
to a suppression of the MET. Nevertheless, the low SM background makes diphoton-l-MET 
a promising search to look for signals of the dark penguin.® 

kinematically allowed, the outgoing photons in Fig. 16 can also be replaced by Z bosons, although the 
smaller BR(vh —^ Xi^) ~ and the additional branching ratios for the Z decays further suppress the cross 
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Figure 16. Dominant process for the diphoton signal. The grey circle indicates the dark penguin. 




Figure 17. Left panel: LHC diphoton bounds for = 2. Right panel: Mt 2 distributions for the 
diphoton signal and background at the luminosity of 3 ab“^, for various values, = ttIxi is 
assumed. For each value of the coupling XVYN is chosen to correspond to a 3cr excess with 
100 fb“^ of data. The perturbative bound on X^/YN depends on the hypercharge of the mediator, 
XVYN <4TrVY. 


The three relevant signal regions for 8 TeV are defined in the first three columns of 
Table 3. For each point in parameter space we compute the constraint from each of the three 
regions, and choose the strongest one. The reach of the 14 TeV LHC is estimated using the 
signal region defined in the last column of Table 3. The results for both 8 and 14 TeV are 
shown in Fig. 17a. 

If a signal is observed, the diphoton+MET channel offers the possibility to measure 
through the stransverse mass variable Mt 2 [81-83]. This variable is constructed from the two 
photons and the missing energy and is defined as 

MT 2 {m) = min [maxima (hn; Pt; Pt )> Pt> Pt )}] > (^-l) 

Pt+Pt=i/>t 

section, making the search for Z hnal states less promising. 
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where 


rriTif^N-, i>T-, Vt ) = Mat + ‘ Pt )• (^-2) 

In the above equations /r at is an unknown trial mass, is the transverse missing momentum, 

is the transverse momentum of photon 7 ^, while the two transverse energies are defined 
as 

$T = = IPt I- (^-3) 

If the trial mass //Af is chosen equal to the xi mass, the distribution of Mt 2 bas an edge at 
the value of In Fig. 17b we show the Mj ’2 distribution at the 14 TeV LHC with 3 ab“^ 
of data, for some illustrative choices of assuming fiN = f^xr ^or each we chose 
the value of the coupling \y/YN such that a 3a excess would be observed at an earlier stage 
of LHC running, namely with 100 fb“^ of data. The edge at Mt 2 ~ can be seen in all 
cases, with some uncertainty of 0(10) GeV due to detector effects. 

In Fig. 17b we assumed the mass of the DM to be known, which allowed us to set 
fiN = ruxi ™ the computation of Mt 2 - But how could be determined experimentally? 
References [84, 85] showed that the value of can be estimated by observing a kink of the 
edge of the Mt 2 distribution as a function of However, a large number of signal events 
is required in order for the kink to be observable, making the application of this method, 
and thus the determination of the DM mass directly in the diphoton process, not likely for 
our model. On the other hand, direct detection constraints may hint that the DM is light. 
As we will discuss in Sec. 7, in the region of parameters where the LHC has sensitivity to a 
dark penguin signal, current direct detection searches require rriy^^ < 10 GeV. Hence one can 
argue that if a diphoton+MET signal is observed at the LHC, and interpreted as involving as 
missing energy the particle that provides the dominant dark matter density in the universe, 
one should use the value /tat ^ 10 GeV. 

It is interesting to compare the results from the monophoton and diphoton searches. The 
comparison requires some further assumption, since the monophoton search sets a constraint 
on the dark penguin coupling between Xh and xu whereas the diphoton search constrains the 
coupling between Xh and Xh- As discussed in Sec. 2 , here we focus on the totally anarchic 
scenario, where all couplings are assumed to be of the same order. In Fig. 18 we show a 
comparison between the 14 TeV monophoton and diphoton searches, for the illustrative choice 
(m^j, = (10, 300) GeV. Diphoton gives a slightly stronger bound, but the difference is 

below the uncertainties associated with our analysis. Notice that, while it is useful to compare 
the two channels within a specific natural scenario such as the anarchic one, in general the 
monophoton and diphoton searches are complementary to each other. 

5.3 Constraint from rare Z decays 

We now turn to the case of light Xh- H the decay Z —)■ XhXh llXiXi is 

kinematically allowed, and is constrained by LEPl data collected at the Z-pole. A search 
performed by OPAL [ 86 ] for events containing at least two photons and missing energy sets 
a strong constraint on this decay. The search region required exactly two photons, satisfying 
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Figure 18. A comparison between the monophoton and diphoton constraints for = 2. The 
perturbative bound on XVYN depends on the hypercharge of the mediator, XVyN < 47rVY. 

the cuts Ej > 1 GeV, \ri^\ < 1.74, > 10 GeV and tt — > 0.0873 (acoplanarity 

angle). In addition, at least one of the two photons was required to have \rj^\ < 1.10, and 
additional photons with > 0.5 GeV were vetoed. No events of this type were observed in 
43 pb~^ of data, corresponding to 1.8 x 10® Z bosons. The expected SM background comes 
from e"*'e“ —?• Z(z^P) 77 , where the two photons arise from initial state radiation, and is quoted 
by OPAL to amount to 0.2 events.^® Therefore, assuming no systematic uncertainty, the 95% 
GL limit is of 2.8 signal events. In order to set a bound on the parameter space, we still need 
to compute the efficiency of the cuts on the signal. For = 30 GeV the efficiency is of 
~ 70%, leading to the 95% CL bound (the branching ratio for the decay Xh —^ Xil is unity, 
since < mz + ) 

BR(Z ^ XhXh) < 2 X 10-® , (5.4) 

which under the assumption of anarchic couplings is much stronger than the other LEP 
constraint from the invisible Z width. The corresponding excluded region is shown in Fig. 19, 
where it is compared to the projected sensitivity of monophoton and diphoton searches at 
the 14 TeV LHC. We conclude that if < m^/2, only the diphoton search will be able 
to marginally improve the LEP constraint. This result motivates our focus on the scenario 
> m^/2, in which case LEP data only give a mild constraint from the invisible Z width. 

^^Simulating the process e''"e“ —>■ in the SM we find a prediction of ~ 0.1 events, in rough agreement 

with the OPAL number. 
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Figure 19. Comparison of the LEP bounds from Z 77 +MET (blue) and from the invisible Z 
width (orange) to the projected LHC sensitivity in monophoton (black) and diphoton (red), at 14 
TeV and with 3 ab“^ of data. The perturbative bound on Xy/YN depends on the hypercharge of the 
mediator, Xy/YN < 47r\/Y. 

6 Non-pointing photon signals 

As discussed in Sec. 2 . 1 , when the flavor structure of the dipole operator is aligned with the 
mass matrix, the decay rate of a heavier fermion Xh into the DM xi a photon can be 
highly suppressed, while the coupling of 7 , Z to a pair of Xh^ or y^’s is still sizable. This 
motivates the study of displaced photon signals from pair produced x^’s, each decaying into 
Xi + 7 - The bounds obtained using the ATLAS non-pointing photon search [77] are shown 
in Fig. 20, for some representative choices of the parameters. The bounds are computed by 
means of a simple counting experiment, by comparing the signal and background yields in two 
different regions, corresponding to |Az-y| > 30 mm (exclusion shaded in blue) and | Az-,,] > 220 
mm (shaded in red). The displacement Az.y was dehned in Fig. 9. For both signal regions, 
an upper bound |Az.y| < 750 mm is also imposed, since the ATLAS paper does not report 
the background expectation for larger values of the displacement. In the first signal region, 
approximately 149 events were expected from the SM background, and 140 were observed. 
In the second region, approximately 14 events were expected and 18 were observed. A few 
comments are in order. The first observation is that, independently of the value of the 
lower bound on the coupling^^ A as a function of e scales as A oc This can be understood 

by noticing from Eq. (3.5) that, for large F^^ (the lower bound on A corresponds to large 
displacement of the decay), the number of signal events passing the cuts is Ns oc A'^F^^, 
since the production cross section scales like A^. Recalling that F^^^ oc A^e^ and that in our 

the following discussion, A is understood as a shorthand for Xy/YN. 
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Figure 20. Bounds in the (e, XVYN) plane obtained from the 8 TeV ATLAS search for non-pointing 
photons. The mass of xi is fixed to 10 GeV, and the mediator mass to 300 GeV. The perturbative 
bound on X\/YN depends on the hypercharge of the mediator, XVYN < 47r\/Y. 


counting experiment the exclusion bound is given by Ng = constant, we obtain the scaling 
A oc Second, comparing the cases = 300 GeV and = 50 GeV we note that 

the bound is weaker in the latter case, especially for the harder cut \^z^\ > 220 mm. The 
reason for this is twofold. In first place, the efficiency of the cuts on the Et of the photons 
and on the (see Sec. 3.5 for details) is smaller for lighter Xh- In second place, a lighter 
Xh is more boosted, which implies that the photon from its decay is typically more collinear 
with Xh- This fact, combined with the upper cut of 4 ns that we impose on the time of 
flight of Xh to avoid spurious collisions due to LHG satellite bunches, implies that the typical 
photon displacement is smaller for lighter Xh- As a consequence, for = 50 GeV the cut 
|Az.y| > 220 mm significantly reduces the signal yield. 

7 Comparison to the direct and indirect detection results 

The comparison between the collider and direct detection constraints is shown in Fig. 21, 
under the assumption of anarchic couplings between the mediators and the dark fermions. 
Here we follow the analysis in [61, 87], comparing the bounds set by the two classes of 
experiments on the magnetic dipole moment defined by the effective operator 

(7.1) 

where x is the DM, and related to the parameters of the dark penguin that couples two DM 
particles to the photon by 

eX'^YN _ 

~ 327r2M ■ 
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Since the momentum exchange in direct detection experiments is small, we can safely integrate 
out the mediators and use the effective description in Eq. (7.1). Results from direct detection 
experiments were obtained in [61], with the assumption of the standard halo model and 
velocities Ugsc = 544km/s, Vq = 232km/s, vq = 220km/s. For the local DM density, the 
standard value p = 0.3 GeV/c^/cm^ was assumed. Figure 21 shows the 90% CL bounds 
from CDMSlite [88], SuperCDMS [89], the XENONIO ^a-only analysis [90], XENONIOO [91], 
LUX [92], CDMS-II-Si [93], CDMS-II-Ge low threshold [93], and CoGeNT2014 data [94], 
together with the 68% and 90% CL allowed regions for DAMA [95] (assuming quenching 
factor Qns = 0.30), CoGeNT2014, and CDMS-II-Si. For XENONIO we take the result with a 
conservative setting of the electron yield to zero below 1.4 keVnr, as in Ref. [91]. For LUX we 
adapt the limit corresponding to zero observed events. For further details on the experimental 
data, as well as on the assumptions made on the low energy thresholds and quenching factors, 
we refer the reader to Ref. [61]. 

We discuss first the = 2 case, shown in the top row of Fig. 21. The red horizontal 
line indicates the upper bound on py. obtained from the invisible Z decay at LEP, whereas 
the green and cyan horizontal lines correspond to the bounds obtained from the projected 
monophoton and diphoton searches at the 14 TeV LHC, respectively. Due to the anarchic 
assumption, all these searches can be interpreted as effectively constraining the size of the 
dark penguin that couples two DM particles to the photon, and thus the size of the effective 
interaction in Eq. (7.1) via Eq. (7.2). This allows us to make the comparison to direct 
detection experiments. 

To gain some understanding of the dependence of the bounds on the mass of Xh: we 
show results using two benchmark values, 50 and 300 GeV. On the other hand, the collider 
bounds are essentially independent of myi in the range considered m-yi < 20 GeV, the only 
appreciable effect being that the LHC monophoton and diphoton bounds worsen slightly for 

> 10 GeV in the case = 50 GeV. The plots in Fig. 21 assume the benchmark value 
of the mediator mass M = 500 GeV, but they serve as useful order of magnitude estimates 
also for heavier mediators. By the end of the high-luminosity LHC run, monophoton and 
diphoton-l-MET searches will be able to test values of the dipole interaction strength that 
are comparable to the current XENONIO reach for rriy^ > 10 GeV, and will provide a much 
better bound than the current direct detection limits for lighter xi- 

The Ny = 1 case is shown in the bottom panel of Fig. 21. The projected 14 TeV 
LHC monojet bound (purple) is slightly worse than the current limit from the invisible Z 
width (red), showing that to improve the constraint it will be necessary to go to future 
colliders. The 100 TeV collider monojet bound is shown in cyan and corresponds to < 
2 X 10“^^ecm, while the ILCIOOO-P monophoton search (green) can reach a bound < 5 x 
10“^®ecm. For comparison, the LEP constraint derived in [42] corresponds to py < 10“^®e 
cm. Similarly to the Ny = 2 case, collider searches play a more important role than direct 
detection experiments when the DM is light, rriy. < 10 GeV. 

The collider bounds can also be compared to the indirect detection searches, where a 
pair of DM particles annihilates into photons. There are two mechanism through which the 
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Figure 21. Top row: comparison between the projected collider constraints and the current bounds 
from direct detection experiments for N-^ = 2, based on the assumption of anarchic couplings between 
the mediators and dark fermions. Here we adapt the direct detection bounds from [61, 87], assuming 
a standard halo model (SHM). See the text for more details. The magnetic dipole moment is defined 
by the effective coupling of the DM to the photon, The red solid lines show the 

bound from the measurement of the invisible Z decay width at LEP. The green (cyan) lines show 
bounds from the monophoton (diphoton) searches at the 14 TeV LHC, with solid (dashed) line for 300 
(3000) fb“^ of data. Bottom: projected collider constraints and current bounds from direct detection 
experiments for iV^ = 1. The lines correspond to the invisible Z width at LEP (red), 14 TeV LHC 
monojet (purple), 100 TeV collider monojet (cyan) and ILCIOOO-P monophoton (green). 


annihilation can proceed in our model: either via two dipole interactions connected by a 
t-channel DM exchange, or via a Rayleigh operator. Using naive dimensional analysis, the 
ratio of the former amplitude to the latter [41] is ~ X'^M/{167r‘^m^). Considering the most 
conservative choice of parameters relevant to our discussion, namely A ~ 4, M ~ 200 GeV 
and ~ 10 GeV, we find the ratio to be about 2. We conclude that for our purposes it is 
a safe approximation to neglect the contribution of the Rayleigh operator, which is typically 
strongly subdominant compared to the diagram with t-channel DM exchange. The latter 
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Figure 22. Upper bounds on the annihilation of DM into photons, derived from the projected 
collider constraints computed in this paper. For = 2 (left) we show the projected LHC bounds, 
while for = 1 (right) the reach of future colliders is presented. In the right panel, the Fermi 
bound corresponds to the expected limit from the 3.7 years data, as reported in [96] (we take the most 
stringent limit in Ref. [96], obtained using the R16 Einasto DM profile). 

dominates the annihilation, with a cross section [40] 

^XX-5'77^ ~ ^ ■ C^-3) 

As shown in Fig. 22, if more than one flavor of dark fermions are accessible at the LHC, the 
size of the annihilation into photons can be tested down to approximately 10“^^ cm^/s. In 
the Ny_ = 1 case, the ILCIOOO-P can set the strongest bound on the annihilation, of order 
10“^^ cm^/s. 

8 Conclusions 

The search for DM plays a central role in the physics program of current and future collider 
experiments. For many plausible DM-SM couplings, a correct description of the DM pro¬ 
duction processes, which is crucial to extract meaningful information from the experimental 
results, requires that the EFT parameterization be UV completed by a simplified model with 
light mediators. This is especially important when the DM-SM coupling arises at loop level. 
In this paper we performed the first simplified model collider study for a loop process, focus¬ 
ing on the Dark Penguin, whose form factors reduce to the magnetic dipole operator at low 
energies. We computed bounds from monojet, monophoton, and diphoton searches at the 8 
and 14 TeV LHC, as well as from the future ILC and 100 TeV hadron collider. Differently 
from searches for EFT interactions, when light mediators are included the optimal search 
strategy requires the cuts on the MET to be as loose as possible, to capture the enhancement 
of the form factors near the threshold for production of on-shell mediators. As we showed 
through a detailed comparison, for light DM mass the collider bounds are complementary to 
those derived from direct and indirect detection experiments. 
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Based on general considerations, it is plausible that the dark sector may be endowed 
with a non-trivial flavor structure. If some of the additional states beyond the DM are 
kinematically accessible, collider experiments offer a unique opportunity to probe the dark 
flavor. By employing a simplified model with a second dark fermion in addition to the DM, 
we showed that, under the natural assumption of anarchic structure in the dark sector, the 
bounds on the DM-SM coupling set by collider searches are much stronger than in the case 
where the DM is the only accessible state. Collider searches not only have the capability to 
probe the flavor structure of the dark sector, but might even allow to measure some of its 
properties. For example, the determination of the mass of the heavy dark fermion could be 
possible by using Mt 2 in the diphoton+MET channel. Furthermore, if the flavor structures of 
the dark penguin and the dark fermion mass matrix are nearly aligned, the decay of the heavy 
dark fermion into DM and photon can be displaced. In this case, it is achievable to extract 
information about the small mixing angle in the dark sector from the search for non-pointing 
photons-l-MET. 

We end with an outlook to future developments. While our study was focused on the dark 
penguin mediating the dipole operator, our method is fully general and can be extended to 
any other loop-mediated interaction between the DM and the SM fields, such as for example 
the Rayleigh and xxG^’^G^i, operators, where is the SM gluon field strength. It would 
also be interesting to analyze in detail the prospects for the ILC, which the preliminary 
results presented here show to be very promising. Additionally, while our first estimates for a 
very high energy hadron collider appear less favorable, the detailed assessment of the design 
requirements that would allow a substantial improvement warrants a dedicated study. 
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A Analytical expressions of the form factors 

In this appendix we present the expressions of the form factors that appear in Eq. (2.2). By 
applying the standard Passarino-Veltman decomposition [97], they can be written as^^ 

_ fq{ui,U2,U3)+ + 

^ [q^ — {mi —mjY][q^ — {mi + mjY]^ ’ 

_ fa{ui,U 2 ,U 3 ) + a^Co{mj,mj,q‘^,M'^,Mj,M‘^) + b^Co{mf, mj, q^, , M^) 

[g^ — {rrii — mj)‘^][q‘^ — {mi + 

(A.l) 

where 

m = Bo{ml M], m2) - Bo{q^ M^), 

U 2 = Bo{ml Ml Ml - Bo{ql M|, M|) , 

U 3 = Bo{m], Mj, Ml - Bo{ml M|, mI . (A.2) 

For the scalar integrals, we followed the convention of LoopTools [98]. The building blocks of 
Fg are given by 

fg{ui,U 2 ,U 3 ) 

= - 2mf{l + U 2 ) - mlq‘^{-4: + ui - 3 u 2 - U 3 ) - A{M‘j - Ml{ui - U 2 - U 3 )] 

- AmiMfq^{ui - U 2 - U 3 ) - 2mj(l + U 2 + U 3 ) - dm?M/(-ui + U 2 + U 3 ) 

- 4m^[-(mi + Mf){ui - U 2 ) + Mjus] - q^[2{M‘j - MI{-ui + U 2 + U 3 ) 

+ q‘^{2 -U 1 +U 2 + U 3 )] - rn‘j[q‘^{-A + ui - 3 u 2 - 3 u 3 ) - 4(M| - MI{ui - U 2 - U 3 ) 

+ 12miMf{—ui + U 2 + U 3 ) + ml—A — 8 ni + 4 u 2 + 6 ^ 3 )] + Amj{mlui — U 2 — 2 U 3 ) 

+ 3m}Mf{ui -U2- U 3 ) + Mfq^{-ui + U 2 + U 3 ) + mj[2M|(ui - U 2 - U 3 ) 

+ 2MI-ui +U 2 + U 3 ) + q^{2u2 + ^ 3 )]}, (A-3) 

agl 2 

= mlmi + Mf) + m‘j[{mi + Mf){2mi + Mf) — 2MI + mlMf{mi + Mf){m’f — 2M|) 

— 2{mi — 2Mf){mi + Mf)Mg — 2M^] — [2mfM/ — 2mjM| + Mj — 2{m^ — miMf + M‘j)Mj 
+ M^]g2 -|- {uii — MlMfq^ + mj{mj — 2Mj + 2Mj — q'^){mi[{mi + Mf){mi + 2Mf) — 2MI 

— Mfq'^} + mj[2mf + AmfMf — miq^ — 2Mf{Mj — Mj + q'^)] + mllmj + AmjMf 

— 2{Mj - M2)2 + ml-2Mj + AMj - 3ql + 2Mjq^ + mi{- 6 Mj + 6 MfMj - AMjql], 

(A.4) 

bg =ag + 2 {mj + mil{—{mj + + 2 Mf)[mjmi{mj + mi) + 2 (m| — mjmi + mlMf] 

+ 4(m2 — mjmi + m-lMj} + 4(mj + mi){2miMf{mi + Mf) + mlmi + 2Mf) — 2miMj 
A-mjlmf + 2miMf + 2{M‘f — Ml]}q^ — 2[mjmi + 2{mj + mi)Mf]q '^. (A.5) 

this appendix, we adopt the shortened notations rriij = m^. . and g"*’® = 
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The building blocks of are 


fa{ui,U2,U3) 

= - 2 m^j - 2 mf + miq^[ 6 {Mj - M^){ui - U 2 - U 3 ) + q ^{-2 + U 1 + U 2 - ua)] 

+ AmjMfq^{ui - U 2 - U 3 ) - + ui + U 2 - U 3 ) + 2 mjmi(-l + U 3 ) 

+ AMfq^{-ui + U2 + U3) + m|[ 2 m ^(2 + U3) - g ^(-4 + ui + U2 + U3)] + m^[ 4 Mjg^(ui 
-U 2 - U 3 ) - 2 mf (-2 + U 3 ) + miq^{A + 3ui + 3^2 + 3 ^ 3 )] + mj{vrifq^{‘i + 3mi + 3m 2 - Sua) 
+ 8miMfq‘^{ui - U2 - U3) - 2mj{l + U3) + g^[6(M| - Mg){ui - U2 - U3) 

+ g^(—2 + ui + U2 + ua)]}, (A-6) 


(X<j 

= — 2{mj — mi)^{mj + mi)^Mg + 2{mj + mi){{mj + Mf){mi + — rrijini + mf 

+ {rrij + mi)Mf — 3Mj] + 2[mj{mi + Mj) + Mf{mi + 3Mf)]Mg — 3M^}g^ — 2[m‘j{mi + 2Mf) 
+ 2Mf{mf + rriiMf — M‘j) — (m* — 2Mj)M^ + mj{mf + 2miMf + 2Mj — M^)]g^ + 2Mjg® , 

(A.7) 


ba = tta — ‘2'{mj — mif‘{mj + mif’ + m* + Mf) — Mg] + 2{mj — mi)‘^{mj + mi) 

X \Mf{mj + mi — 2Mf) + 2Mf] + 2{mj + mi) [Mf{mj + ruj + 3Mf) — 3Mf] — 2Mjg® . 

(A.8) 


The above expressions can be easily implemented in a numerical code, and provide stable 
results for all values of g^, in particular for g^ > 4M^ (where M denotes generically the 
mediator mass), where the form factors develop an imaginary part as a consequence of the 
virtual states going on-shell. For completeness, we also report the expressions of the form 
factors in terms of integrals over Feynman parameters, which are more suitable for analytical 
expansions. These read 


where 



(A.9) 


Af = — xyq^ — zxmf — yzm'j + {x + y)M‘j + zMg , 

A, = AfiMf^Ms), 

Cq = y{2y + z- 2)mj - {y + z - l){2y + z)mi - {2y + z - l)Mf , 

dq = y{2y + z - l)mj - {y F z - l){2y F z - l)mi + {2y + z - l)Mj , 

Co- = - yzmj + z{y + z-l)mi +{z-l)Mf , 

da = — yzmj + z{y + z — l)mi — zMf , (A.10) 
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and we defined f d?x = dxdydz 5{x + y + z — 1). It is straightforward to check numerically 
that the Eqs. (A.9) agree exactly with Eqs. (A.l) for < 4M^, whereas for > 4M^ the 
Feynman parameter integrals can become numerically unstable, and only the expression in 
terms of scalar integrals should be used.^^ 

By using Eqs. (A.9), one can derive the expressions of the form factors for large mediator 
mass, Mf = Mg = M ^ rrii, mj, 


1 

’ 



(A.ll) 


In this limit, the amplitude in Eq. (2.2) can be seen as generated by the effective Lagrangian 


I 

1 + dij 


{X^hg'YN ^ 

I 927 r 2 M 2 Bf,, + h.c. 


+ 


1 {X%g'YN_ 


1 + 6ii 647r2M 


+ h.c.. (A.12) 


B Statistics 


In this appendix we briefly describe our procedure for setting limits on the DM parameter 
space. Given the expected number of background events Nb, we exclude a signal model 
yielding Ns events if 

/ Nobs sp 

— {x;NB + Ns)dx<p, (B.l) 

where dP/dx {x-,Nb + Ng) is the normalized probability distribution function for signal plus 
background, N^bs is the number of events observed, and p is the chosen probability. For 
example, for a 95% CL exclusion, p = 0.05. Notice that we are setting a ‘one-sided’ limit. 
Throughout our analysis we neglect systematic uncertainties on the signal, because they are 
subleading to those associated with the background. 

Under the assumption that signal plus background follows a Poisson distribution with 
mean Nb + Ng, and neglecting all systematic uncertainties, the exclusion limit is given by 


r(A^ofes + 1; 0) 


(B.2) 


where r(s,( 7 ) = P~^e~^dt. For example, for p = 0.05, we find for Nb = N^bs = 0 that 
Ng ~ 3.0. 

If we assume that signal plus background follows a Gaussian distribution with mean 
Nb + Ng and standard deviation Ng + (dNB)"^, where 6 Nb = \JNb + {eNB^, with e the 
relative systematic uncertainty on the background, we find the exclusion limit 


^ ^erf-^ (1 - 2p). (B.3) 

have cross-checked our results against Ref. [49], where the form factors were computed for the special 
case rrii = rrij, Y = 1/2, N = 2. We find agreement, except for a few small differences that we believe are due 
to typos in App. A of [49]: 1) in their Eq. (A7), a factor z should multiply the second term on the right-hand 
side (RHS), and an overall factor (—1) should multiply the RHS; 2) in their Eq. (A8), an overall factor (—2) 
should multiply the RHS. 
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For example, for p = 0.05 the right-hand side is equal to ~ 1.645. 

Equation (B.2) was used to set the limit from prompt diphoton and from the rare Z 

decay Z —>■ 77 -l-MET, whereas all the remaining limits in this paper were computed using 

Eq. (B.3). 
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